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[erical  evaluation  of  the  e- integral  occurring  i7t 
the  theolorsen  arbitrary  airfoil  potential  theory 

Irven  Naiman 

SUMMA3  " 

A  more  precise  ra     I  of  evaluating  the  e-integral 
occurring  in  the  arbitrary  airfoil  theory  :f  Theodorsen 
(  t\CA   heps.  No  a .  411  and  452)  is  developed  by  retaining 
/    her  order  terraa  in  the  Taylor  expanaion  and  by  use 
of  Simpson's  rule.   Formulas  are  riven  for  routine  calcu- 
lation of  the  e -integral  and  for  the  necessary  compu- 
tational coeff icients.   The  computational  coefficients 
are  tabulated  for  a  40-point  division  of  the  range  of 
integration  from  0  t">  2tt.   With  no  increase  in  compu- 
tational work  the  systematic  error  in  the  numerical 
value  of   e   is  reduced  from  the  order  of  1  percent  to 
approximately  0.1  percent. 

I]  TRODUCTI  ' 


The  solution  of  the  general  problem  by  means  of 
conformal  transformation  for  the  flow  about  an  arbitrary 
airfoil  (references  1  and  2),  a  syi  etrica]  lattice 
(reference  3),  and  a  biplane  (     cence  4)  involves  t]  t 

determination  of  the  imaginary  part  of  a  cor-]  lex  trans- 
formation function,  given  the  real  part.   As  ~' !  wn  in 
references  1  and  kJ  the  real  part  may  be  e     ded  in  a 
Fourier  series  and  the  imaginary  part  is  the  conjugate 
Fourier  series.    It  is  else  shown  in  these  references 
that  the  imaginary  part   e  may  be  obtained  from  t]  e 
real  part   \!/  by  the  following  functional  equation: 

r2TT 
e((pi)  =  _^-  /   *(<p)  cot  *   "  CP'  do       (1) 
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EVALUATION  0?  TIC  e- INTEGRAL 


The  evaluation  of  the  ^-integral  is  complicated  by 
the  discontinuity  at   cr  =  o'  .   This  difficulty  nay  be 
surmounted  by  a  separate  solution  across  the  discon- 
tinuity.  When   s  =  cp  -  cp'   is  substituted  in  equation  (1) 


£(<?')  =  -57? 


i     r2rf-w  3 

/      h{^?  +  s)  cot  ; 


•cp' 


or,  because  of  the  periodicity  of  this  function, 


e(cp')  = 


..2-  / 


2tt- 


'|r(cp'  +  s1)  cot  2    ds 


(2) 


The  discontinuity  now  occur?  at   s  =  0.   Tor  purposes  of 
numerical  evaluation  this  integral  ma^  be  broken  up  as 

f  cllov  ?  r 


e(cpt)  =  -i-  J    \J/(cp'  +  s) 

-J  -s 


cct  ~- ds  +  i     \|/(p'  +  s)  cot  -1  d: 


J. 


ei  +  £? 


(3) 
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v.'her? 


rz 


■1   - 


1       I 
2 


i  \ 

i 

-'-a 


«>(cp?  +  s)  cot  %   ds 


(4) 


md 


r2TT-a 


"  1 


l|ffi>«  +  s)  cot  £  ds 


(5) 


Evaluation  of   c'i  •  ~  The  first  integral   c,   includes 

the  discontinuity  and  the  limit   s  may  be  taken  as  some 
convenient  small  value.   By  a  Taylor  °eries  expansion  the 
integral  is  easily  evaluated  as  follows: 


i!'!v'  +  a)  =  l'(v')  +  s\|/'  (q>«)  +  4-i|r»(<pi)  +  |^|fiii(<pi)  +  ... 


When  this  expansion  is  substituted  in  equation  (4)  the 
internals  containing  the  even-ordered  derivatives  are 
found  to  be  identically  zero.   Equation  (4)  then  becomes 


_1_ 
2rr 


.'3 


C  0  Xj    7t    \jl  s 


\|/"»         i 


7.  I 


CO1 


ds   + 


(6) 


where  the  derivatives  are  evaluated  at   cp' .   The  Taylor 

g 

expansion  for   cot  -^   is 


cot  -^ 


3 
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ind  equation  (6)  is  thus  obtained  as 


fi  = 


;rr 


■R'-"  \-  36  "3600  ~  •••/ 


Vs 


/ 


s°\l/"»  [  1 


£0   1680 


o5    V  /       Co  2 

150 T  V    84 


'J      8820*    V^   108   '  "  -y   ' 


7  ; 


Evaluation  of   £o.-  The  second  integral   Co  (equa- 

tion  (5))  nay  be  rearranged  for  convenience  in  numerical 
calculation  as  follows: 


"2tt 


~TT 


/    \|/  (cp '    +   s)    cot  75  ds  +    I 

-'s  ♦% 


,-2rr-s 


^r(cp'  +  s)    cot  §  ds 


_1_ 

"2tt 


i  r 


'S 


•vMcp'  +  s)    cot  ^   ds+     /    \|/(c?'  -  s)    cot  -  ds 


m 


"1 

v|/(co'  +  s)    -  \J/(p»  -  s)  j  cot  |  ds  I 

J 


( ■■■■) 


where   -s   has  been  substituted  for   s   in  the  second 
integral  and  the  Units  have  been  rearranged  accordingly. 


numerical  ksth< 


[  of  reference  1„~  In  reference  l  the  interval 

0  to   2tt   is  divided  into  n   equal  parts  of  magnitude 

2rr/n   (n   is  an  even  number)  .   The  values  of   \|/  are 
d<  signated  \j/  «.   \!r  „   ,..<,,  \J/ 


Vn,       \ 


J-       *1>        '  '  '  > 


w]  ere     ^0     it    the   value   of     \j/     at      cp  =   cp<      and     V:   n  =  U.-n 

"2  2 


o.  L4D27 


the  value  at  cp  =  c?»  ±  rr.   The  integrations  are  per- 
ked over  intervals  of  width  Srr/n  with  the  ty     values 
the  midpoint  of  the  interval.   The  range  of  integra- 
tion for   e-i   i3  frc;n   s  =  -Tr/n   to   s  =  rr/n  and  for 
e0   from  njn      to   2ir  -  — . 

The  first  integral   c^   i-s  evaluated  by  retaining 
only  the  first-order  terns  in   s, 


el  —^V  =  -|V  (;) 


where  bhe  slope  vv '   is  determined  graphically  a1 
that  is,  at   cp  =  <p ■  (s  =  0). 

The  second  integral   £o   is  composed  of  the  si 

integrals  across  each  interval 


2k+l, 

n-1  r~--- 

£o  -  - 

c"   f~:  i'2k-l 

k=l  •• TT 

n 


^    '<tt 


1  ^     n 
7—    /       I       *(<Pf  +  s)  cot  £  ds 


The  function  v!/  does  not  change  much  across  the 
interval  and  is  therefore  approximated  by  use  of  the 

cp '  +  — —  ).   Then , 
n  / 


n-1     c — ——it 


n 


n  1  •   2k+l 

,  n-1  sin  —  ; — rr 

-  -  -J  ^>   M/v  log ^_ 

7~i  sm  — '- — TT 

K=l  2n 
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or,   by  equation    (8),    with     ^n_^ 


\b  i 


n 


/2 


NT 


k 


k=l 


(10) 


where 


ak  =    log 


n 


=    0 


..      2k  +   l 

cm  — - Tr 

_2ri 

2k   ••    1 

s  m  — — — — rr 
2n 


(II) 


?he    complete   integral   is    given  by      €   ~    ci  +    €2>      or 

£  -  -<§V  +f/L   ak(^  "  *-*)f      (12) 

L  ~k=1  J 
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rhe   Newton-Stirling  formulas  for  derivatives    (refei 
ence    5,    p.    75)    are 


1     *,  1         E 


gS^fll 


C.  -T 


-   53'"    - 


**  * 


■'   -  '    =    6%    - 


J 


(13) 


ore   s  =  tabular  In4       (2-Tr/n)  .   The  mean  central  dif- 
ferences  6'lc  can  be  expj  of  the  tabular 

values   as 


26^      =    \|/n     -    ty 


-1 


o,  I ,         ,      o»  I 


■-1    "      -: 


(14) 


=   \I/-z    -    4tyo    + 


2   +   sxi;^   -   5V_i   +   4V_2    -      V_3 


J 


gives 


substitution  of  relations  (13)  in  eq   I  *  sn  (7) 


£i  =  - 


4, 

L 


g^«6'H2 


s 

13 


. 


-  25 


'^  ~ 


O.J  *  "3 


^  /  13 


1512 


o  •   /   "   • 


(15) 
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The    further   substitution   of   relations    (14)    in   equa- 
tion   (15)    gives 


el   = 


J-  Ai, 
"2i\\cl 


'-l) 


(, 


19 


2    +  ■£   +  ■—■  + 


90     •  •  -y 


4/x  +  iL  +  A 

18\  15        G4 


\^c        V-2;ilg  +oo5  +  ...1  -  S     ^&70  +  37q   +  •  .  .^    - 


+  (^3-^3) 


19 


450  +  *•  7  "  S"l 


f  1 

x^j. . 


+ 


'J 


i 
J 


or 


el   ~   blO-l    ~  *l)    +   b2(V-2    ~   Uf2)  +   t3(-i/_3   +   Mr5)     (16) 


where 


bl   =  ffl 


/        4      19\  p/l  i  5     \ 

y      9      90/   "    s    ^18      135      1512/ 


1  /239         167        2 

'    2ttV90          2520    3 

-    ... 
/ 

^E-^r(Mi)- 

2/11 

S    U?0        378 

1   /83            2        2 

Tn\225    ~    315   s 

\ 

h            1  /l9             s2 

/ 

/ 

°c    =    2ir\450        1512 

1 


>    (17) 
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e    3ecor.d   integral      r-0      is    evaluated  by  Simpson's 

rule    from     \i/-,       to      ^n._i 


c2   -    -:•-'—  ~-!'-'i     cot   ~    +    4\!/c    cot 


2ir 


.t   or,  I    1 


11 


+    :'v: ,    cot   —   - 
n 


or,    by  cquatS  •/.•,,•: 


n/ 


/o 


wttere 


CO 


TT 


c  .    kir 

■5—   cot   -  - 
Oil 


cn    E    0 


^ 


(19) 


where    (except   in   the    first   term)      c  =  2      for     k     odd  and 
a  -   4      for     k      even.        The    complete   integral    is   gj 

c  =    Cl   +    £2,      or 


2Z? 


j..  -   X 


(20) 


10 
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where 


At,     -     1>-     +     Ci 


!1) 


Values  of   Av   for  n  =  40   are  Riven  in  table  I. 


ACCURACY  0?  EVALUATION 


The  accuracy  of  the  two  methods  of  evaluation 
described  may  be  determined  by  integrating  various 
harmonics.   The  results  are  presented  as  ratios  of  the 
integrated  value  to  the  correct  value  so  that  a  valve 
of  unity  is  a  correct  evaluation.   Values  of  this  ratio 
for  the  harmonics  are: 


Harmonic   j   40-point  method 

of  reference  1 

4  0-point  pre  sent 
method 

1 

2 
3 
4 
5 
10 

1.01434 
1.02662 
1 .03700 

1.04  547 
] .05280 

1.03112 

0.99977 

I  -CO ')44 
r   :  04.0 

1.00060 
.9  J  8  iO 

.95~::0 

Inasmuch  as  the  higher  harmo:  ics  • ~     ~. 
in  a  much  smaller  proportion  than  the 
3uch  contours  as  are  encountered 
error  01  i ;' ■  ".•  40-point  method  of 
order  of  1,    rcent,  whereas  tr. 
presence'  -   ?in  is  approximately  0.1  D.eroe: 


'  :j.nter 


for 


■ s -  the 


is  y. 


the 


met] 
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TABLE   I.-    VALUES    OF      Ak      FOR  USE  WITH  EQUATION    (20) 


k 

A 

1 
:  ^40   j 

J 

H: 

n  =  4  ) 

1 

1 
0.62827 

!  n 

0,01423 

o 

.14824 

12 

.02422 

5 

.  v,'614 

1? 

.01021 

4 

. 10259 

14 

.01698 

5 

.04024 

!   15 

.00690 

6 

.06542 

1   16 

.■'1083 

i 

.02720 

17 

.00400 

8 

.04E  i 

18 

•  OOP  2 8 

9 

.0.1951 

1   ] 

.00066 

10 

. 03333 

:   20 

'  .            

0 

i 

TABLE   II.-   VAU~r    DP      ak      PC  :   USE  WITH   EQUATION    (12); 

;:::r":  or  i  ^        ss  l  and  2 


ak 

'  '  1 

k 

n  =  10 

n  =  20 

n  =  40 

1 

1.06544 

1.090::'" 

1.09056 

2 

.44211 

.49426 

.50671 

3 

.23117 

.31141 

.33028 

4 

.10302 

.21750 

.24303 

5 

0 

.1577? 

.19028 

6 

=11449 

.1545? 

7 

.0! 

.12  041 

8 

.05115 

.10326 

9 

.02493 

.03207 

10 

0 

.07362 

11 

.06714 

12 

.05711 

13 

.048ie 

14 

.04004 

15 

.  ^,  \J  i-J  \.>   I  1 

16 

.  Kj  <0  O  *_/ 1_/ 

17 

.01837 

18 

. 01245 

19 

.00613 

20 
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